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Abstract 

We study experimentally nodal domains of wave functions (electric field distributions) lying in the 
regime of Shnirelman ergodicity in the chaotic microwave half-circular ray-splitting rough billiard. 
For this aim the wave functions of the billiard were measured up to the level number N = 415. 
We show that in the regime of Shnirelman ergodicity (N > 208) wave functions of the chaotic 
half-circular microwave ray-splitting rough billiard are extended over the whole energy surface 
and the amplitude distributions are Gaussian. For such ergodic wave functions the dependence 
of the number of nodal domains on the level number N was found. We show that in the 
limit N — > oo the least squares fit of the experimental data yields ^n/N ~ 0.063 ± 0.023 that is 
close to the theoretical prediction ^n/N ~ 0.062. We demonstrate that for higher level numbers 
N ~ 215 — 415 the variance of the mean number of nodal domains cr%/N is scattered around the 
theoretical limit cr^/N ~ 0.05. We also found that the distribution of the areas s of nodal domains 
has power behavior n s oc s~ T , where the scaling exponent is equal to r = 2.14 ± 0.12. This result 
is in a good agreement with the prediction of percolation theory. 
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In recent theoretical papers by Bogomolny and Schmit [1] and Blum et al. 2] the distribu- 
tions of the nodal domains of real wave functions $f(x, y) in 2D quantum systems (billiards) 
have been considered. Nodal domains are regions where a wave function ty(x,y) has a def- 
inite sign. The condition ty(x,y) = determines a set of nodal lines which separate nodal 
domains. Bogomolny and Schmit [1] have proposed a very fruitful, percolationlike, model 
for description of properties of the nodal domains of generic chaotic system. Using this 
model they have shown that the distribution of nodal domains of quantum wave functions 
of chaotic systems is universal. Blum et al. [2J have shown that the systems with integrable 
and chaotic underlying classical dynamics can be distinguished by different distributions of 
the number of nodal domains. In this way they provided a new criterion of quantum chaos, 
which is not directly related to spectral statistics. 

n n 

Theoretical findings of Bogomolny and Schmit [1] and Blum et al. [2| have been recently 
tested in the experiment with the microwave half-circular rough billiard by Savytskyy et al. 

In this paper we present the first experimental investigation of nodal domains of wave 
functions of the chaotic microwave ray-splitting rough billiard. Ray-splitting systems are a 
new class of chaotic systems in which the underlying classical mechanics is non-Newtonian 
and non-deterministic . In ray-splitting systems a wave which encounters a disconti- 

nuity in the propagation medium splits into two or more rays travelling usually away from 
the discontinuity. Ray splitting occurs in many fields of physics, whenever the wave length is 



large in comparison with the range over which the potential changes. Idea 



the investigation of ray-splitting phenomena are ray-splitting billiards 



cavities with dielectric inserts 5|, 8,1^, llOj. Measurements of wave functions of ray-splitting 



model systems for 
3] and microwave 



systems are very demanding because in principle they require the direct access to the all 
parts of the system ljj including those filled with ray-splitting media, such as dielectric 
in the case of ray-splitting microwave billiards. This is one of the main reasons for which 
only low wave functions (N < 100) of ray-splitting billiards have been measured so far 111 ]. 
In this paper we use a new method of the reconstruction of wave functions introduced by 
Savytskyy and Sirko [3] which in the case of the half-circular microwave ray-splitting rough 
billiard allowed for the reconstruction of wave functions with the level numbers N < 415. 

In the experiment we used the thin (height h = 8 mm) aluminium cavity in the shape of 
a rough half-circle (Fig. [1]) which consisted a half-circular Teflon insert of radius Rd = 8.465 
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FIG. 1: Sketch of the chaotic half-circular microwave ray-splitting rough billiard which consists 
a half-circular Teflon insert of radius Rd = 8.465 cm. Dimensions are given in cm. The cavity 
sidewalks are marked by 1 and 2 (see text). Squared wave functions \%>n(R c , 9)\ 2 were evaluated 
on a half-circle of fixed radius R c = 19.25 cm. Billiard's rough boundary is marked by T. 



cm. The insert had the same height as the rough cavity. The microwave cavity simulates 
the rough ray-splitting quantum billiard due to the equivalence between the Schrodinger 
equation and the Helmholtz equation 6]. This equivalence remains valid for frequencies 
less than the cut-off frequency v c = c/2r]h ~ 13.1 GHz, where c is the speed of light and 
t] = 1.425 is the index of refraction of the Teflon insert. 

The cavity sidewalls were made of two segments. The rough segment 1 is described by the 
radius function R(9) = R + Y,m=2 a m sm(m9 + 4> m ), where the mean radius i?o=20.0 cm, 
M = 20, a m and m are uniformly distributed on [0.084,0.091] cm and [0,27r], respectively, 
and < 9 < n. It is important to note that we used a rough half-circular cavity instead of a 
rou gh c ircular cavity because in this way we avoided nearly degenerate low-level eigenvalues 
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141 ] . Additionally, a half-circular geometry of the cavity was necessary for the accurate 



measurements of the electric field distributions inside the billiard. 

According to [15] the roughness of a billiard may be characterized by the function k(9) = 
(dR/d6)/ Rq. The roughness parameter k defined as the angle average of the function k(9) 
was for our billiard k = {(k 2 {9)) g) 1 ^ 2 — 0.200. In such a billiard the dynamics is diffusive 
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in orbital momentum due to collisions with the rough boundary because the roughness 
parameter k is much larger the chaos border parameter k c = M -5 / 2 = 0.00056 [l^l- The 



roughness parameter k determines also other properties of the billiard 16j. The eigenstates 
are localized for the level number N < N e = 1/128A; 4 = 5. The border of Breit-Wigner 
regime is given by N w = M 2 /A8k 2 ~ 208. It means that between N e < N < N w Wigner 
ergodicity [3] ought to be observed and for N > Nw Shnirelman ergodicity should emerge. 
In the regime of Shnirelman ergodicity wave functions have to be uniformly spread out in 
the billiard 17j. In this paper we focus our attention on Shnirelman ergodicity regime. 



It is worth noting that rough billiards and related systems are of considerable interest 
elsewhere, e.g. in the context of microdisc lase rs |l8 . ljj, light scattering in optical fibers 20]. 
ballistic electron transport in microstructures jsnj, dynamic localization 2^] and localization 
in discontinuous quantum systems [231 ] . 

In order to measure the wave functions (electric field distributions inside the microwave 
billiard), which are indispensable in investigation of nodal domains, we used a new, very 
effective method described in |12| . It is based on the perturbation technique and construction 
of the "trial functions" . 

Following 12J we will show that the wave functions ^/^(r,6) (electric field distribution 



En(t,9) inside the cavity) of the billiard can be determined from the form of electric field 
E N (R C , 6) evaluated on a half-circle of fixed radius R c (see Fig. [I]). 



The first step in evaluation of E^{R C , 6 , 
the perturbation technique developed in 



is measurement of \En(R c , 0) \ 2 - For this purpose 



24j and used successfully in 



24, m 



26 



27] was 



applied. In this method a small perturber is introduced inside the cavity to alter its resonant 
frequency according to 

v -v N = v N {aB 2 N -bE 2 N ), (1) 

where is the iVth resonant frequency of the unperturbed cavity, a and b are geometrical 
factors. Equation (1) can be used to evaluate Ej^ only when the term containing magnetic 
field Bn is sufficiently small. In order to minimize the influence of on the frequency shift 
v — vn a small piece of a metallic pin (3.0 mm in length and 0.25 mm in diameter) was used as 
a perturber. The perturber was attached to the micro filament line hidden in the groove (0.4 
mm wide, 1.0 mm deep) made in the cavity's bottom wall along the half-circle R c and moved 
by the stepper motor. Application of such a small pin perturber reduced the largest positive 
frequency shifts to the uncertainty of frequency shift measurements (15 kHz). It was verified 
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FIG. 2: Panel (a): Squared wave function \^4i5(R c ,d)\ 2 (in arbitrary units) measured on a half- 
circle with radius R c = 19.25 cm (^415 — 12.98 GHz). Panel (b): The "trial wave function" 
^415(^,0) (in arbitrary units) with the correctly assigned signs, which was used in the recon- 
struction of the wave function ^415 (r,6) of the billiard (see Fig. [3]). 



that the presence of the narrow groove in the bottom wall of the cavity caused only very 
small changes 5un of the eigenfrequencies i/jy of the cavity \5vn\/vn < 10~ 4 . Therefore, 
its influence into the structure of the cavity's wave functions was also negligible. A big 
advantage of using the perturber that was attached to the line, was connected with the fact 
that the perturber was always vertically positioned, which is crucial in the measurements of 
the square of electric field E^. The influence of the thermal expansion of the Teflon insert 
and the aluminium cavity into its resonant frequencies was eliminated by stabilization of 
the temperature of the cavity with the accuracy of 0.05°. 

The regime of Shnirelman ergodicity for the experimental rough billiard is defined for iV > 
208. Using a field perturbation technique we measured squared wave functions \^ n(R c , 9)\ 2 
for 30 modes within the region 215 < N < 415. The range of corresponding eigenfrequencies 
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was from 1/215 — 9.42 GHz to ^415 ~ 12.98 GHz. The measurements were performed at 0.36 
mm steps along a half-circle with fixed radius R c = 19.25 cm. This step was small enough to 
reveal in details the space structure of high-lying levels. In Fig. [2] (a) we show the example 
of the squared wave function \^n(R c , 9)\ 2 evaluated for the level number N = 415. The 
perturbation method used in our measurements allows us to extract information about the 
wave function amplitude \^ n{Rc,§}\ at any given point of the cavity but it doesn't allow to 



determine the sign of ^n{Rc,@) 29] . However, the determination of the sign of the wave 



function ^/ N (R C , 6) is crucial in the procedure of the reconstruction of the full wave function 

nn 

\I/jv(r, 6) of the billiard. The papers [3|, |12f] suggest the following sign- assignment strategy. 
First one should identify of all close to zero minima of \^n(R c , 8)\. Then the sign "minus" 
is arbitrarily assigned to the region between the first and the second minimum, "plus" to 
the region between the second minimum and the third one and so on. In this way the "trial 
wave function" tyisf(R c ,6) is constructed. If the assignment of the signs is correct the wave 
function \I/jv(r, 9) should be reconstructed inside the billiard with the boundary condition 
^v(r r ,#r) = 0. 

The wave function of a rough ray-splitting half-circular billiard outside of the half-circular 
Teflon insert (r > R d ) may be expanded in terms of Hankel functions 

L 

*T{r, 0) = £ a s n s (k N r) sin(s0), (2) 

8=1 

where fi a (ar) = Re{H® (x) + S ss (k N R d )H^(x)} and k N = 2-kv n /c. H^(x) and Hf\x) 
are Hankel functions of the first and the second kind, respectively. The matrix S ss '(k^R d ) 



is defined as follows 



30] 



HW(k N R d ) - V [Ji( V k N R d )/J s ( V k N R d )]HP(k N R 



HP (k N R d ) - ^J'XvknR^/JsivkNRdW^ikNRd) 
where the derivatives of Hankel and Bessel functions are marked by primes. In Eq. (2) 
the number of basis functions is limited to L — kj^r max + 3, where r max = 20.7 cm is the 
maximum radius of the cavity. l™ ax = k^r max is a semiclassical estimate for the maximum 
possible angular momentum for a given k^. The functions with angular momentum s > l^ ax 
describe evanescent waves. We checked that the basis of L wave functions was large enough 
to properly reconstruct billiard's wave functions. The coefficients a s may be determined 
from the "trial wave functions" ^n(R c ,8) via 

a s = £n s (k N R c )}- 1 f V N (R C , 0)sin{s9)d0. (4) 
l Jo 
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FIG. 3: The reconstructed wave function ^415 (r,9) of the chaotic half-circular microwave rough 
billiard. The amplitudes have been converted into a grey scale with white corresponding to large 
positive and black corresponding to large negative values, respectively. Dimensions of the billiard 
are given in cm. The position of the half-circular Teflon insert of radius Rd = 8.465 cm is marked 
with a solid line. 

The wave functions of the billiard inside the Teflon insert (r < R d ) may be expanded in 
terms of circular waves 

v 

*N(r, 0) = £ a' s J s (r]k N r) sin( S 0). (5) 

s=l 

In Eq. (5) the number of basis functions was limited to V = rjk^Rd. The coefficients a s 
given by Eq. (4) and the continuity condition fulfilled at the border of the dielectric insert 
^N^iRdiO) = ^ l N{Rd,9) may be used to evaluate the coefficients a' s in Eq. (5) allowing in 
this way to reconstruct the full wave function ^(r, 6) of the billiard. 

In the evaluation of the coefficients a' s in Eq. (5) an important role plays the value of the 
refraction index n of the Teflon insert. We measured the refraction index 77 = 1.425±0.002 of 
Teflon by measuring the set of resonant frequencies of a microwave circular cavity of radius 
Rt = 3.25 cm entirely filled by it. 

Using the method of the "trial wave function" we were able to reconstruct 30 experimental 
wave functions of the rough half-circular billiard with the level number N between 215 
and 415. The wave functions were reconstructed on points of a square grid of side 4.2 • 

7 



10 -4 m. As the quantitative measure of the sign assignment quality we chose the integral 
7 / r \iffjtf(r, 9)\ 2 dl calculated along the billiard's rough boundary T, where 7 is length of 
T. In Fig. [2] (b) we show the "trial wave function" \E f 415 (i? c , 9) with the correctly assigned 
signs, which was used in the reconstruction of the wave function ^415 (r, 9) of the billiard (see 
Fig. [3]). It is worth noting that inside of the Teflon insert the size of nodal domains are much 
smaller than outside of it. The remaining wave functions from the range N = 215 — 415 
were not reconstructed because of the accidental near-degeneration of the neighboring states 
or due to the problems with the measurements of \^n(R c , 9)\ 2 along a half-circle coinciding 
for its significant part with one or several of the nodal lines of \I/7v(r, 9). The problem of the 
near-degenerated states is important because in the presence of the perturber the resonances 
are shifted, which may cause the initially non-overlapping states to become near-degenerated 
at certain positions of the perturber. Such a situation prevents us from the reconstruction 
of the wave functions. The problems mentioned are getting much more severe for iV > 200. 
Furthermore, the computation time t r required for reconstruction of the "trial wave function" 
scales like t r oc 2™ z_2 , where n z is the number of identified zeros in the measured function 
\^n(R c ,9)\. 



The structure of the energy surface 15| of the billiard's wave functions plays an important 
role in the identification of their ergodicity. To check it we extracted wave function ampli- 
tudes C^P = (n, l\N) in the basis n, I of a half-circular ray-splitting billiard (desymmetrized 



annular ray-splitting billiard) 28j with radius r max and a half-circular Teflon insert of radius 



Rd . The normalized eigenfunctions of the half-circular ray-splitting billiard are given by 

x , M Aj n J/(»7Kj n r)sin(Z0), < r < R d , 

®ni(r,9) = < (6) 

[ A in [Ci n Ji(n ln r) + Di n Yi(Ki n r)\ sin(/0), R d < r < r maX} 

where A ln = {f (/<f d rJ^K^rfdr + r[G n J z (/w) + D ln Yi{K ln r)f dr))~* . Ji(K ln r) 
and Yi(Ki n r) are Bessel and Neumann functions, respectively. The main quantum number 
n = 1, 2, 3 . . . enumerates the zeros yi n = Ki n r ma x of the radial function 

CmMym) + D ln Yi(yi n ) = 0, (8) 

and / = 1, 2, 3 ... is the angular momentum quantum number. The coefficients C/ n and Di n 
can be determined from the continuity conditions of the wave function $ n ;(r, 9) and it's 
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FIG. 4: Structure of the energy surface in the regime of Shnirelman ergodicity. Here we show 
the moduli of amplitudes |C^j | for the wave functions: (a) N = 215, (b) N = 415. The wave 
functions are delocalized in the n, I basis. Full lines show the energy surface (see text). 



The moduli of amplitudes | and their projections into the energy surface for the rep- 
resentative experimental wave functions N = 215 and N = 415 are shown in Fig. HI As 
expected, in the regime of Shnirelman ergodicity the wave functions are extended over the 
whole energy surface [l3j|. The full lines on the projection planes in Fig. Hl(a) and Fig. H^b) 
mark the energy surface of a half-circular annular ray-splitting billiard H(n, I) ~ En = k N 
estimated from the formula \H(n,l) — En\/En < 0.12. The peaks |Cj^| are spread almost 
perfectly along the lines marking the energy surface. 



derivative & nl (r, 8) on Teflon's boundary R d 






FIG. 5: Panel (a): The amplitude distribution P^nA 1 / 2 ) for the wave function N = 215. Panel 
(b): The distribution P^nA 1 / 2 ) for the wave function N = 415. The amplitude distributions 
were constructed as histograms with bin equal to 0.2. The width of the distribution P{^>) was 
rescaled to unity by multiplying normalized to unity wave function by the factor A 1 / 2 , where A 
denotes billiard's area. Full lines show standard normalized Gaussian prediction Pq(^ A 1 / 2 ) = 
(l/v^e-* 2 ^ 2 . 

Ergodic behavior of the measured wave functions can be also tested by evaluation of 



the amplitude distribution P(^f N ) [3l|, |32j. For irregular, chaotic states the probability of 
finding the value at any point inside the billiard should be distributed as a Gaussian, 
P(^n) ~ e - ^*^. In Fig. H3(a) we show the amplitude distribution P(ty ^A 1 ^ 2 ) for the 
wave function N = 215 while in Fig. 0(b) the distribution P{^> ^A 1 / 2 ) for the wave function 
N = 415 is presented. The distributions were constructed as normalized to unity histograms 
with the bin equal to 0.2. The width of the amplitude distributions P($n) was rescaled to 
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FIG. 6: The number of nodal domains H^r (full circles) for the chaotic half-circular microwave ray- 
splitting rough billiard. Full line shows the least squares fit Hat = a\N + b\^/N to the experimental 
data (see text), where ai = 0.063 ± 0.023, b\ = 0.77 ± 0.40. The prediction of the theory of 



Bogomolny and Schmit 



l|] ai = 0.062. 



unity by multiplying normalized to unity wave functions by the factor A 1 / 2 , where A denotes 
billiard's area (see formula (23) in 31]). For all measured wave functions lying in the regime 
of Shnirelman ergodicity the distributions of P^nA 1 ^ 2 ) were in good agreement with the 
standard normalized Gaussian prediction Pq(^>A 1 / 2 ) — (1/ V / 27r)e _,I ' 2A//2 . 

The number of nodal domains K^r vs. the level number N in the chaotic microwave ray- 
splitting rough billiard is plotted in Fig. [61 The full line in Fig. [6] shows the least squares 
fit K N = aiN + hy/N of the experimental data, where a x = 0.063 ± 0.023, b x = 0.77 ± 0.40. 
The coefficient a\ = 0.063 ± 0.023 coincides with the prediction of the percolation model 
of Bogomolny and Schmit [1] N^v/iV ~ 0.062 within the error limits. The errors of the 
coefficients a± and b\ are relatively high because the number of nodal domains fluctuates 
significantly in the function of the level number N, what was also demonstrated in [2] (see 
Fig .(5)). It is worth mention that in the paper [3J the coefficient a\ was estimated in the 
experiment with the microwave rough billiard without the ray-splitting Teflon insert. Its 
value a\ = 0.058 ± 0.006 was also close to the theoretical prediction. The second term 
in the least squares fit corresponds to a contribution of boundary domains, i.e. domains 
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FIG. 7: The variance of the mean number of nodal domains divided by the level number <J%/N for 
the chaotic half-circular microwave ray-splitting rough billiard. Full line shows predicted by the 
theory limit cr^/N ~ 0.05, Bogomolny and Schmit 



that include the billiard boundary. Numerical calculations of Blum et al. Q performed for 
the Sinai and stadium billiards showed that the number of boundary domains scales as the 
number of the boundary intersections, that is as \^N. Present results together with the 
results of {3} clearly suggest that in the rough billiards (with and without ray-splitting), at 
low level number N, the boundary domains also significantly influence the scaling of the 
number of nodal domains Kjv, leading to the departure from the predicted scaling ~ N. 

Measured wave functions of the ray-splitting billiard may be also used for the calculations 
of the variance a 2 N of the mean number of nodal domains. It was predicted in [1] that for 
chaotic wave functions the variance of the mean number of nodal domains should converge 
to the theoretical limit <j 2 n ~ 0.05iV. In Fig. [7] the variance of the mean number of nodal 
domains divided by the level number u%/N is shown for the microwave ray-splitting rough 
billiard. The variance a% = N 1 _ 1 J2^=i(^Ni ~ (^n)) 2 wa s calculated in the window of N w = 5 
consecutive eigenstates measured between 215 < N < 415, where the mean number of nodal 
domains was defined as (Kzv) = j^- Y.i=i ^at 4 - For level numbers N < 300 the variance cr%/N 
is above the predicted theoretical limit, however, for 300 < N < 415 it is slightly below it. 
A similar erratic behavior of <J 2 N /N was also observed in [l|. 
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FIG. 8: Distribution of nodal domain areas. Full line shows the prediction of percolation theory 
log 10 ((rt s /n)) = -^f logio((s/s m j n }). The least squares fit log 10 ((n s /n)) = a 2 - r \og 10 ((s / s m in)) 
of the experimental results lying within the vertical lines yields the scaling exponent r = 2.14±0.12 
and a% = —0.06 ± 0.12. The result of the fit is shown by the dashed line. 



The percolation model [lj allows for applying the results of percolation theory to the 
description of nodal domains of chaotic billiards. The percolation theory predicts that the 
distribution of the areas s of nodal clusters should obey the scaling behavior: n s oc s~ T . The 



scaling exponent 33j is found to be r = 187/91. In Fig. [8] we present in logarithmic scales 
nodal domain areas distribution (n s /n) vs. (s/s m i n ) obtained for the microwave ray-splitting 
rough billiard. The distribution (n s /n) was constructed as normalized to unity histogram 
with the bin equal to 1. The areas s of nodal domains were calculated by summing up 
the areas of the nearest neighboring grid sites having the same sign of the wave function. 
In Fig. [8] the vertical axis (n s /n) = ■j^ : J2i f =i n i N ^/ n ^ N " > represents the number of nodal 
domains n^ N ' of size s divided by the total number of domains averaged over Nt = 30 
wave functions measured in the range 215<iV<415. In these calculations we took 
into account only the nodal domains which entirely lied outside or inside of the Teflon 



insert for which percolation theory 33j should be applicable. The horizontal axis in Fig. [8] 
is expressed in the units of the smallest possible area s^ n [lj], (s/s m i n ) = ^J2^i s / S mim 
where s^2 n = ^(joi/t^aO 2 an d joi — 2.4048 is the first zero of the Bessel function Jo(joi) = 0. 
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For nodal domains lying inside the Teflon insert the refraction index was according to our 
measurements r\ = 1.425 while outside of the insert we assumed 77 = 1. The full line in Fig. [8] 
shows the prediction of percolation theory log 10 ((n s /n)) = — ^ log 10 ((s/s mi „)). In a broad 
range of log 10 ((s/s m j n )), approximately from 0.2 to 1.4, which is marked by the two vertical 
lines the experimental results follow closely the theoretical prediction. The least squares fit 
log 10 ((n s /n)) = a 2 — rlog 10 ((s / s min )) of the experimental results lying within the vertical 
lines gives the scaling exponent r = 2.14 ± 0.12 and a 2 = —0.06 ± 0.12, which is in a good 
agreement with the predicted r = 187/91 ~ 2.05. The result of the fit is shown in Fig. [8] by 
the dashed line. 

In summary, for the first time we measured high-lying wave functions of the chaotic 
microwave ray-splitting rough billiard. We showed that in the limit N — > 00 the least squares 
fit of the experimental data reveals the asymptotic number of nodal domains Kjv/iV ~ 0.063± 
0.023 that is close to the theoretical prediction K^/iV ~ 0.062 We demonstrate that for 
higher level numbers N ~ 215—415 the variance of the mean number of nodal domains a 2 /N 
is scattered around the theoretical limit a 2 /N ~ 0.05. Following the results of percolationlike 
model proposed by [l] we confirmed that the distribution of the areas s of nodal domains 
has power behavior n s oc s~ T , where scaling exponent is equal to r = 2.14 ± 0.12. This 
result is in a good agreement with the prediction of percolation theory [33j, which predicts 
t = 187/91 ~ 2.05. The experimental results presented in this paper strongly suggest that 
many properties of nodal domains in chaotic ray-splitting billiards are the same, like in 
conventional chaotic billiards without ray-splitting. 
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